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The magnetoresistance of thin superconducting (SC) strips subject to a perpendicular magnetic 
field B and low temperatures T manifests a sequence of alternating SC-insulator transitions (SIT). 
We study this phenomenon within a quasi one-dimensional (ID) model for the quantum dynamics 
of vortices in a line-junction between coupled parallel SC wires, at parameters close to their SIT. 
Mapping the vortex system to ID Fermions at a chemical potential dictated by B, we find that a 
quantum phase transition of the Ising type occurs at critical values of the vortex filling, from a SC 
phase near integer filling to an insulator near 1/2-filling. For T — >• 0, the resulting magnetoresistance 
R(B) exhibits oscillations similar to the experimental observation. 
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The conduction properties of low-dimensional super- 
conducting (SC) systems (thin films and wires) are 
strongly dominated by fluctuations in the SC order pa- 
rameter. A particularly prominent manifestation of the 
role of fluctuations is the appearance of a finite dissi- 
pative resistance below the mean-field critical tempera- 
ture T c of the bulk superconductor. At low temperatures 
T <C T c , the dominant fluctuations are in the phase of 
the complex order parameter. Most notably, topological 
excitations (vortices and phase-slips) can generate dis- 
sipation in their liquid state. In the T — > limit, their 
quantum dynamics becomes significant and may drive a 
transition to a metallic or insulating state 0, Q • 

In the one-dimensional (ID) case, i.e. SC wires of 
width and thickness smaller than the coherence length £, 
the resistance essentially never vanishes at finite T due 
to thermal activation of phase-slips 0, B (for T < T c ) 
or their quantum tunneling at lower T 0, [a, • In con- 
trast, in the 2D case (SC films) superconductivity is well- 
established at sufficiently low T. However, a quantum 
(T —> 0) superconductor-insulator transition (SIT) [l], Q 
can be tuned by an external parameter which leads to 
proliferation of free vortices. Employing charge-flux du- 
ality 8] it is possible to view the SC phase as a vortex 
solid, and the insulator as a vortex superfluid. 

A convenient means of inducing a SIT in SC films 
is by application of a perpendicular magnetic field B. 
At fixed T, a positive magnetoresistance R{B) is typi- 
cally observed in a wide range of B. The SIT is then 
clearly indicated in the data as a crossing point of these 
isotherms at a critical field B c , separating a SC phase 
(where dR/dT > 0) for B < B c from the insulating phase 
(dR/dT < 0) for B > B c . 

A recent experimental study of a strip geometry [9( 
- namely, a SC wire of width comparable to £ - offers 
an opportunity to probe the crossover from a ID to 2D 
quantum dynamics of the topological phase-defects in SC 
devices. The prominent observation is that in the pres- 
ence of a perpendicular field B, the magnetoresistance 
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FIG. 1: (color online) Isotherms of R as a function of 
B, for J = 1.06K, V = 0.98K, V-/L = 0.01K [see Eq. 
(|12[) ]. indicating superconducting (SC) and insulating (I) sec- 
tions alternately. Inset: R vs. the gap A<j [see Eq. (|15p 
and text therafter] near a single critical point, for T = 
0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5K. 



R(B) exhibits oscillations which amplitude is sharply in- 
creasing at low T, in striking resemblance to the behavior 
of Josephson arrays [l(| and SC network systems [Tlj . 
Moreover, the SIT at B c appears to be preempted by 
several consecutive transitions at lower fields, from a SC 
to an insulator or vice versa alternately. 

The periodicity of the above mentioned oscillations is 
consistent with a single flux penetration to the sample, 
suggesting that the observed SC or insulating behavior 
of the system is determined by commensuration of vor- 
tices within the strip area. When an integer number of 
vortices can be fitted along the strip length, superconduc- 
tivity may be supported even at sufficiently high B such 
that a large fraction of the sample area turns normal. De- 
viation from commensurability of the vortex filling weak- 
ens superconductivity, possibly inducing a transition to 
a metallic [l(| or insulating state. 

In this paper we focus on the strongly quantum fluctu- 
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ation regime characterizing the ID vortex matter close to 
an SIT, and propose a theory for its low T transport be- 
havior. The system is shown to exhibits a series of quan- 
tum phase transitions of the Ising type, manifested as 
SC-insulator oscillations of the Ohmic resistance R(T, B) 
[see Fig. [Tj. This result underlines a correspondence 
between charge-flux duality across a SIT and the order- 
disorder duality characteristic of the Ising transition. 

We consider a SC strip subject to a perpendicular mag- 
netic field B < B c . Assuming that the high vortex den- 
sity in this case leads to near merging of their cores along 
the central axis of the strip, we model the system as a 
line-junction formed by a pair of parallel SC wires of 
length L 3> £, separated by a thin normal barrier of 
width w. In the low T (phase-fluctuations) regime, the 
dynamics of the collective phase field in the wires (4>i (x, t) 
with i = 1, 2) is governed by the effective ID Hamiltonian 



Hq — Hi + Hi + Hi, 



in which (using units where K = 1) 



Hi = 



dx 



UoPt 



4m 



(1) 



(2) 



Hi, 



dx [~gj cos(0i - (f> 2 - qx) + U p\p 2 ] ■ (3) 



Here the operator pi{x) denotes density fluctuations of 
Cooper pairs in wire i, and can be represented as (l2j ] 

Pi {x) = --dJiix) +PoY e * 2 P(WM*)) ( 4 ) 

7T 

in terms of the conjugate field Oi(x) satisfying 
[<j} t {x),d x 6 z (x')} = inS(x' - x). The first term in Eq. © 
hence describes a charging energy; p s is the superfluid 
density (per unit length) assumed to be monotonically 
suppressed by increasing B, po = Ps(B — 0) and m is 
the electron mass. The inter-wire coupling [Eq. ([3])] 
consists of a Josephson term and an inter-wire Coulomb 
interaction, of coupling strengths gj and U, respectively; 
q = 2ttwB/$q (with <I>o the flux quantum) denotes vor- 
tex density per unit length. Hq describes an ideal system, 
to which we later add a disorder scattering potential. 

It is convenient to introduce total and relative phase 
fields via the canonical transformation cj)± — {(pi ± 
4>2)/^/2 and the corresponding conjugate fields ~d x 6±, 
in terms of which the Hamiltonian ([T]) is separable: 



H Q =H + + H_ 



(5) 



where H+ = H^ 



H_=H 



(-) 

LL 



dx 



-gj cos(-\/2</)_ — qx) + g c cos(\/8#-) 



H 



(±) 

LL 



V± 

2tt 



dx 



K±(d x e ± )+— (dM 2 

K± 



(0) 



and the parameters are given by 
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(7) 



Here we have accounted for the most relevant interaction 
terms, neglecting umklapp terms that are effectively sup- 
pressed due to the rapidly oscillating factor in Eq. (|4|. 
We next define new canonical fields 



-^=0_ , 9 = V29_ 



(8) 



in terms of which H LL acquires the form of a Luttinger 
Hamiltonian with an effective Luttinger parameter K = 
K-/2. Assuming that the original parameter K- is close 
to the critical value for a SIT in a ID wire, K c = 2 0, 
we obtained K 1. This yields 



H_ 



V- 

2^ 



[{d x ef + (d^) 2 ] (9) 

+ / dx [—gj cos(20 — qx) + g c cos(26>)] . 



The model Eq. (fTU)) can be refermionized by introduc- 
ing right (R) and left (L) moving Fermion fields [l3T ] 



_ e ±ik F x e i(T<l>+0) 



(10) 



in terms of which H_ becomes a free Hamiltonian. Note 
that the "Fermi momentum" kp is dictated by q, which 
can be viewed as a vortex filling factor in the line-junction 
area. Quite interestingly, this implies that near a SIT, 
it is natural to adapt a duel representation of this sys- 
tem in terms of fermionic vortex fields. We note that 
the vortex-matter in the system is constrain by an effec- 
tive periodic potential set by the vortex- vortex interac- 
tion and the boundary conditions: the leads connected 
to x — ±L/2 (assumed to be macroscopic superconduc- 
tors), and the strip edges which induce an effective "im- 
age charges" potential [3]. Hence, the vortices tend to 
form a pinned chain where core positions are separated 
by a constant spacing a = L/N, in which N — T[qL] 
(I[z] the integer value of z) denotes the total number of 
vortices ■ We hereon regard Eq. (fT0"|) as a continuum 
limit of a lattice model, where the coordinate x = na (n 
integer). Consequently, q is replaced by the deviation of 
the vortex density from the closest commensurate value: 

w{B-B N ) $ 
q = , B N = NB , B = — ; 11) 

N = I[B /Bq] so that < q < |. The lattice constant a 
also determines the short-distance cutoff a in Eq. (fT0"|) . 
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The fermionic representation of H_ is given by 



has the same form but with modified velocities and gaps. 
In particular, the modified gaps for finite fi v are given by 



= J dx{v-[rp R (x)(-id x )rp R (x)-ip f L (x)(-id x )'tp L (x)] A u4 (fj, v ) 



J±V , V= y/V 2 + 



(15) 



- fi v [ip R (x)ip R (x) + ipl(x)ip L (x)] 

- J[^ R {x)i) L {x) + 4>l(x)yj R (x)} 

+ V[^ R {x)^ L {x) + (12) 

where J — nagj, V — nag c and the vortex chemical 
potential fi v is determined by the band-structure of the 
Fermions. For simplicity, we associate the commensu- 
rate case q = with a 1 /2— filled band at fi v = 0; e.g., 
fi v = to cos[7r(l/2 — q)} with e ~ v~/a. We note that the 
particular dependence of [i v on q is not crucial to demon- 
strate the qualitative features of the model: however gen- 
erally, it is monotonically increasing with fi v (q = 0) = 0. 

Following the analogous problem of spin- 1/2 ladders 
(l3l HH , it is useful to decompose the complex Fermions 
[Eq. (|10j) ] in terms of the Majorana fields 
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iV2 



ty u - 41) (13) 



(y = R,L). Recasting Eq. (TT21 in fc-space and using the 



Fourier transformed fields 



f]ju,k — Vjv,-k 



here A^ = J ± V denote the gaps in the excitation 
spectrum for commensurate vortex filling (q = fi v =0), 
in which case %k decouples into two independent blocks. 
Since J, V are positive, the u sector is higher in energy. 

We now focus on the case of interest, where the system 
is assumed to be in the SC phase but close to a SIT so 
that the Josephson energy J is slightly larger than V, and 
-C . In this case, the high energy sector u can 
be truncated, and the low-energy properties are governed 
by the <i-type Fermions. In particular, the gap can 
change sign upon tuning of J below the critical value 
J c = V where A^ = 0. Indeed, for fx v — each species 
of free massive Fermion models described by (Q3]) can be 
independently mapped to an Ising chain in a transverse 
field @,[l6j]. When finite vortex "doping" is introduced 
by tuning B away from Bn such that fi v ^ 0, the original 
d and u sectors mix. However, the resulting long wave- 
length (fc — > 0) theory can still be cast in terms of two 
decoupled d and u sectors, where the energy spectrum 
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, we obtain 
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-v-k ) 



V 

( VlR,k , V2L,k , r]2R,k , T]lL,k ) 5 (14) 



While A„ remains positive and large for arbitrary fi v , a 
quantum phase transition occurs at a critical value of ji v 
[which can be traced back to a sequence of critical fields 

B^ NS> via fJ- v (q) and Eq. (fTTj) ]. where A^ changes sign. As 
B -t B ( C N \ | A d | ~ \B-B ( C N) \. Below we show that these 
Ising-like quantum critical points correspond to SIT. 

To study transport properties of the system it is neces- 
sary to include a scattering potential, generically induced 
by random, uncorrelated impurities along the coupled 
wires. Without loss of generality, we hence consider ran- 
dom impurities in wire 1 by including a linear coupling of 
Pi(x) to a disorder potential Vd{x) in the Hamiltonian. 
The leading contribution to dissipation arises from the 
backscattering term of the form 



H D = J dx((x) cos{26i(x)} 



(16) 



[see Eq. (TJ|], where we assume (((x)) — 0, (((x)((x')) — 
DS(x — x') [with () including disorder averaging]. 

For sufficiently small D, a perturbative treatment of 
Hd (see, e.g., Chap. 7.2 in [lij]) yields the d.c. resistance 
at finite T 



/>oo 

R(T) kAD dtt%m{ X (t)} 
Jo 



x (t) = (sin{20i(t)} sin{20 1 (O)}) o ; (17) 

A is a numerical constant, and ()o denotes an expec- 
tation value with respect to H [Eq. (|5J]. Using 9\ = 
-^(0 + + 9-), the correlation function x(t) decouples into 

x (t) = xc+(t)xc-(t) + xs+(t)xs-(t); 

xc ± = (cos{V29 ± (t)}cos{V29 ± (0)}) ± (18) 

XS ± = (sm{V29±(t)}sm{V26 ± (0)}) ± 

where ()± are evaluated w.r.t. H±. Since H+ is a Lut- 
tinger Hamiltonian [see Eq. (j^J], we obtain 12 1 



Xc+(t) = Xs+(t) = lim 



sinh{7rT (t - ie)} 



(19) 



In contrast, as discussed below, Xc-(t)> Xs-(t) depend 
crucially on the parameters of (|14p. and in particular on 
the magnitude and sign of the masses A Ui< j. 

To evaluate xc-i Xs- we first note that in terms of the 
field 6 [Eq. ([8])]. they correspond to correlation functions 
of cos (9, sin#, which lack a local representation in terms 
of Fermion fields. However, a convenient expression is 
available in terms of the two species of order (cr U)C j) and 
disorder (a u ,d) Ising fields 13j,ll6|: for A^ > 0, 



cos* 



(20) 
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For Ad < 0, the roles of ad, (id are simply interchanged. 
The correlators xc-j Xs- carL therefore be expressed in 
terms of C x (t) = (a x (t)a x (0)), C x (t) = (a x (t)a x (0)) 
(A = u, d), which have known analytic app roximations 
in the semi-classical regime (|A A | » T) OEEll: 

C x {t) ~ lA^KoiilA^t), C x (t)^\A x \^ 4 (21) 

[with Ko(z) the modified Bcssel function]. In the quan- 
tum critical regime (| A d | < T), C d (t) ~ C d (t) - t" 1 / 4 . 

Employing Eqs. ffTT ]) (f2T]l . we derive expressions for 
the low-T resistance R(T, B) near commensurate fields 
Bn [Eq. (|TT|) ] where = A^, and near -Bat+i 
A + i)£>o, where A^ is maximally negative. Neglecting 
terms of order e~ A "/ T , we obtain for B ~ Sjv 



i?(T,B) a i? s 



A d {B) ( TraAd(B) 



-k(B) 



-A d (B)/T 



T \ 



= 2-KT 1 (B), R s ozD{A d (B)y 



(22) 



Superimposed on a moderate increase with _B arising 
from k(B) due to the suppression of p s [Eq. ([7J)], the 
exponential factor leads to a strong decrease and i? — > 
at T — > as long as A^(i?) > is finite. The disordered 
Ising phase is thus identified as superconducting, suggest- 
ing that the fields <7d physically represent phase-slips. In 
contrast, for B ~ B N+ i (Ad < 0) we find 



fi (T,B)^(-l^)i) 1/4 (^)' (23) 

(i?» cx D). Since if" 1 < |, n(B) > yielding i? -> oo 
at low T, indicative of an insulating behavior. In fact, 
in this regime the perturbative calculation leading to Eq. 
([23| is not valid in the T — > limit, where localization 
takes over and R(T) diverges exponentially 12j. 

The above analysis implies that the quantum critical 
points at (where A^ = 0) correspond to SC-I and 

I-SC transitions alternately (see Fig. 1). However, note 
that unlike the 2D SIT, R(T) in their vicinity does not 
manifest a metallic behavior: the power-law correlations 
in the critical regime yield 

R(T, B) ~ !T^ S ) , 1 {B) = i - k(B) < . (24) 

This reflects a slightly moderated insulating behavior, an 
asymmetry that stems from the quasi-lD nature of the 
system. As a result, a sharply defined crossing point of 
isotherms does not exist, as clearly indicated in Fig. 1. 

To summarize, we have shown that prominent fea- 
tures of the magnetoresistance oscillations in SC strip- 
like devices are captured by a toy-model for the quantum 
Josephson-vortices in a line-junction between parallel SC 
wires. When the wires are close to a SIT, this system can 
be described by a field theory of ID free Fermions, im- 
plying the existence of a sequence of quantum critical 



point of the Ising-type. It should be pointed out that 
since the Fermions are massive, small deviations from 
our ideal choice of parameters such that the Fermions 
become interacting do not change the essential proper- 
ties, as the interactions can be treated perturbatively. 
The low-T behavior of the resistance R(T, B) provides a 
transparent interpretation of the critical points as transi- 
tions from a SC (near integer vortex fillings) to insulator 
(near ^-integer vortex fillings) or vice versa. However, 
its behavior in the vicinity of the transitions is distinct 
from the SIT in fully 2D SC films. In particular, it does 
not reflect the duality symmetry of the underlying model, 
and the insulating regimes are effectively widened. 

Final note added: During the preparation of the 
manuscript we became aware of an independent work 
[l9| . proposing an alternative theoretical model as in- 
terpretation to the data of Ref. 0. Both models yield 
qualitatively similar magnetoresistance oscillations. 
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